Discovery of X-ray Emission from the Wolf-Rayet star WR 142 of 

oxygen subtype 



ON 
O 

o 

(N 
£ 

00 
(N 

P? 

Or 
6 

o3 



> 

m 
in 
in 

o 

ON 
O 



13 



L. M. Oskinova, W.-R. Hamann, A. Feldmeier 
Institute for Physics and Astronomy, University Potsdam, 14476 Potsdam, Germany 
lidaOastro . physik . uni-potsdam . de 

R. Ignace 

Department of Physics and Astronomy, East Tennessee State University, Johnson City, TN 31614, USA 

Y.-H. Chu 

Department of Astronomy, University of Illinois, 1002 West Green Street, Urbana, IL 61801, USA 



ABSTRACT 

We report the discovery of weak yet hard X-ray emission from the Wolf-Rayet (WR) star 
WR142 with the XMM-Newton X-ray telescope. Being of spectral subtype W02, WR142 is a 
massive star in a very advanced evolutionary stage, short before its explosion as a supernova or 
7-ray burst. This is the first detection of X-ray emission from a WO-type star. We rule out 
any serendipitous X-ray sources within ps 1" of WR142. WR142 has an X-ray luminosity of 
Lx = 7 x 10 30 ergs _1 , which constitutes only <10 -8 of its bolometric luminosity. The hard X- 
ray spectrum suggests a plasma temperature of about 100 MK. Commonly, X-ray emission from 
stellar winds is attributed to embedded shocks due to the intrinsic instability of the radiation 
driving. From qualitative considerations we conclude that this mechanism cannot account for the 
hardness of the observed radiation. There are no hints for a binary companion. Therefore the 
only remaining, albeit speculative explanation must refer to magnetic activity. Possibly related, 
WR 142 seems to rotate extremely fast, as indicated by the unusually round profiles of its optical 
emission lines. Our detection implies that the wind of WR 142 must be relatively transparent to 
X-rays, which can be due to strong wind ionization, wind clumping, or non-spherical geometry 
from rapid rotation. 
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Stars of the Wolf-Rayet (WR) type have a 
highly peculiar chemical composition and very 
strong stellar winds. The WR spectra are sorted 
into three subclasses: WN, WC, and WO, ac- 
cording to the dominance of nitrogen, carbon, 
and oxygen emission lines, respectively. These 
spectral types cor respond to evolutionary stages 



type stars, that represent the final evolutionary 
stage of a massive star prior to its explosion as 



a typ e Ic supernova or 7-ray burst (jHirschi et al 
l2005h . 



(jConti et al.1 11983). The largest C+O abundance 



and the fastest stellar winds are observed in WO 



The WR stars continue to c hallenge our un- 
derstanding of line-driven winds. Schaerer ( 1996?) 
pointed out the importance of the iron opacity for 
the acceleration of WR winds. The first hydrody- 
namical model for a WR w ind was presented by 
iGrafener fc Hamannl (|2005f ) . In their simulation 
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the mass loss is initiated at high optical depth by 
the so-called "iron bump" in the opacity. It was 
thus demonstrated that WR-type winds can be 
driven by radiation pressure. 

It has long been known that line-driven winds 
are subject to an instability that ca n lead to 



strong shocks ( Lucv fc White 1980h . These 



shocks are thought to explain the X-ray emis- 
sions from O star winds, as predi cted by time- 



dependent hydrodynamic m odeling (jOwocki et al 



1988tlFeldmeier et al.lll997l) and largely confirmed 
by observations (Kramer et al.ll20 03: Osk inova et al 



2006 | IZhekov fc Pallall2007t IWaldron fc Cassinellil 
20071 ) . The gro wth of instability in WR w inds was 
investigated by iGavlev fc Owocki ( 1995 ). They 
found that despite of damping effects due to the 
multi-line scattering, the instability remains effec- 
tive. Therefore, X-ray emission from wind shocks 
could, in principle, be expected in WR winds, a 
conjecture that has not been yet tested by time- 
dependent hydrodynamic simulations. 

Significant observational effort has b een made 

to stu dy the X-ra y emissi on of WR stars. IWhite fc Lone 



1986) 



Pollocld (|!987h . iPollock et al 



(1995). 



Oskinoval (|2005l ) presented X-ray observations 
of Galactic WR stars. A survey of X-ray emis- 
sion from WR star s in the Magella n ic Clou ds 
was conduc t ed by Guerrero fc Ctru (l2008al rbh. 
Ignace et akl (|200ot ) and lOskinoval (|2005f ) demon- 
strated that X-ray properties of single WR stars 
differ from those of O stars. Whereas O stars dis- 
play a trend in which the ratio of the X-ray to the 
bolometri c luminosity Lx /Lhn] h as a typical value 



of 1(T 7 (TLong fc White I Il98d iBerghoefer etaL 
1997 ; Sana et akl 2006f ). this trend is not observed 
in the case of WR stars. 

Observations with the XMM-Newton and Chan- 
dra X-ray telescopes established th at some bona 



fide sin g le WN stars are X-r ay active ([Skinner et al 



single . 

2002al|rJ llgnace et al.l l2 o"o3lOskinovall2005l). while 
other s are apparently not (IOskinoval2005l : Gosset et al 



20051 ). lOskinova et al.l (|2003f ) found that no sin 



gle WC star had been conclusively detected at 
X-ray energies, a result that continues to hold. 
Among the WR subclasses, only the WO-type 
stars have not been observed in X-rays so far. In 
this Letter we present the XMM-Newton observa- 
tions of the closest WO type star WR142. The 
star is introduced in Sect. 2, and its XMM-Newton 
observations are described in Sect. 3. The impli- 



cations of the XMM-Newton detection of WR 142 
are discussed in Sect. 4. 

2. The WO-type Star WR142 

WR142, also named Sand 5 and St 3, has a 
spectr um characteristic for the spectral subtype 



WQ2 (jBarlow fc Hummedll982tlKingsburgh et al 
19941 ). T he optical spec t rum o f WR142 was dis- 



cussed bv IPolcaro et al. ( 1997 ). who also noticed 
some line variability. 

Figurc[T] shows the spectral energy distribution 
(SED) of WR 142 together with model calculated 
with the Potsda m Wolf-Rayet (PoW R) model at- 
mosphere code ( Grafener et al.lr2002h . Photomet- 
ric IR measurements are plotted together with 
our optical spectrum and the Spitzer IRS mid- 
IR spectrum. We have adopted parameters typ- 
ical for a WO star: stellar tem perature = 
150 k K, "transformed radius" (cf. iGrafener et al. 
20021 ) R t = 2R & , and a composition of 40% 



carbon, 30% oxygen and 30% helium (by mass). 
WR 142 is assumed to be a member of the open 
cluster Berkeley 87 at a distance of d = 1.23kpc 



( Turner et al. 20061 ). Based on this preliminary 



model, the fit of the photometric observations 
requires an interstellar reddening of Eb-v — 
1.7mag and a stellar luminosity of logLboi/^o = 
5.35, implying a stellar radius of only i?„ = 
0.6 i?©. The corresponding mass-loss rate is about 
10 _51 AfQyr^ 1 for an adopted microclumping vol- 
ume filling factor of 0.1. 

The adopted model does not provide an en- 
tirely satisfactory fit to the line spectra. For ex- 
ample, the model does not match the huge Ovi 
emission at 3811, 3834 A , a pr oblem also experi- 
enced bv lCrowther et all ( 2000f ) while reproducing 
these lines for Sand 2 with the CMFGEN code. We 
also fitted the SED shown in Fig. Q] with a model 
that has a mass-loss rate lower by a factor of two, 
10~ 5 ' 4 M Q yr~ 1 , and higher bolometric luminosity, 
logLboi/-^© = 5.65. This model fits the SED in 
Fig.Q] equally well. Figure[2] shows the radius of 
unity optical depth plotted as function of wave- 
length in the X-ray range for both models. The 
stronger wind is opaque even to hard X-rays, but 
the thinner wind is largely transparent, because 
its higher ionization reduces the X-ray absorbing 
ions. The same could happen in denser, but hotter 
models. Thus our ability to predict the influence 
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of wind photo-absorption is somewhat limited, ow- 
ing to ambiguities in the ionization state of metals 
and uncertainty in the mass-loss rate. 

The profiles of the emission lines in the spec- 
trum of WR 142 are very broad. Assuming that 
the line widths correspond to the wind termi- 
nal velocity, the velocity of w 5500 km s -1 
would be deduced. However, the profile shapes 
of almost all lines are much more round than the 
roughly Gaussian shapes usually seen in WR spec- 
tra. Such round profiles cannot be reproduced 
by the standard models. It is tempting to repro- 
duce these profiles by convolution with the semi- 
ellipse for rotational broadening, albeit rotating 
stellar winds certainly require a more sophisti- 
cated treatment which has not been accomplished 
yet. If rotation is the cause for the round pro- 
files, the projected rotation speed must be com- 
parable to Voo, i.e. the star would rotate near 
to break-up with v TOt sinz<4000kms _1 (see inset 
in Fig. [I]). Interestingly, a similar suggestion has 
been made for the hottest and the most compact 
Galactic WN star, WR2, (T* = 140kK,i^ = 



0.89 R e?,, Vrnt sini w 2000km s ) (jrlamann et al 



2006) 



3. Observations 

WR 142 was observed by XMM-Newton during 
two consecutive satellite orbits (Obsld 0550220101 
and Obsld 0550220201). The data were merged 
and analyzed using the latest versions of software 
SAS 8.0.0. After the high background level time 
intervals have been rejected, the combined expo- 
sure time of all detectors was « 100 ks; an EPIC 
image of WR 142 and its surroundings is shown 
in Fig.[3J WR142 is detected with a 5a confi- 
dence level in all XMM-Newton EPIC detectors 
using standard source-detection algorithms. We 
define a "hardness ratio FIR" as HR = (N^ard — 
N S oft)(N h ard + Ngoft)" 1 , where N soft is the number 
of counts in the 0.25-2.0 keV band and Nhard m the 
2.0-12.0 keV band. For WR 142 we find HR« 0.57. 
For comparison, the hardness ratio of the WN star 
WR 1, which has a reddening similar to WR 142, is 
HR(WRl)w -0.9, while more reddened WR110 
has HR(WR110)« -0.4 and less reddened WR6 
has HR(WR6)« —0.6. Thus X-ray emission from 
WR 142 is the hardest among putatively single 
WR stars. The count rates in the 0.25 keV- 



12.0keV band are (1.89 ± 0.34) x lO^cs" 1 for 
EPIC MOS1+2 and (3.80 ± 0.84) x lO^cs" 1 for 
EPIC PN cameras. Assuming a two-temperature 
thermal plasma model (fcTi = 0.3 keV, kT^ = 
lOkeV), the observed X-ray flux of WR142 is 
F x = 4 ± 2 x lO-^crgs^ciTT 2 . The redden- 
ing towards WR 142 is known from the analysis 
of its optical spectrum, and the distance is known 
from its cluster membership. The X-ray luminos- 
ity of WR142 is thus L x « 7 x lO^ergs^ 1 , or 
logLx/^boi « -8. 

The angular resolution of XMM-Newton is < 6". 
To exclude the potential confusion with a source 
in close vicinity of WR142, we inspected opti- 
cal and infra-red images with higher angular res- 
olution. Accordin g to the USNO-B1.0 catalog 
(|Monet et al.ll2003h . the closest object to WR142 
is located 8" away. The optical monitor (OM) on 
board of XMM-Newton provides images with an 
angular resolution of rj 1". The OM images in 
four filters did not detect any objects closer than 
8" around WR142. The Spitzer space telescope 
has observed the Berkeley 87 cluster in the infra- 
red (IR). With an angular resolution of < 1", the 
Spitzer IRAC camera took images in four chan- 
nels. No point sources within 8" from WR 142 
were detected in any of the IRAC channels. It is 
extremely unlikely that any X-ray source is not de- 
tectable in either the optical or IR. We therefore 
rule out any serendipitous X-ray sources within 
« 1" of WR 142 that might get confused as our 
target. 

In addition, there is no evidence to suggest 
that the observed X-rays originate from a wind 
blown bubble around WR 142. Although theoreti- 
cally expected, there is a dearth of WR stars with 
detected diffuse X-ray emiss i on from their wind 



blow n bubbles (|Chu et alJ 120031 : IWrigge et al 



2005). The only two detected hot bubbles show 



a limb-brightened morphology and are extended 
on the scale of parsecs. On the contrary, we con- 
fidently detect a point source at the position of 
WR142. 

The spectra of WR 142 were extracted from the 
15" region. The small number of counts obtained 
in our observation does not allow for quantitative 
spectral analysis, but interesting conclusions can 
be made about the gross energy distribution. Fig- 
ured] shows the X-ray spectrum of WR 142 before 
a background subtraction, together with the spec- 
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trum of a nearby background region normalized to 
the same area. The emission from WR 142 domi- 
nates over the background at 2 — 7 A (1.7 — 6 keV) , 
while longwards it dives below the background un- 
til about w 20 A (0.6 keV) where it rises above the 
background again. The presence of "a dip" be- 
tween w 7 — 12 A (0.8 — 1 keV) is also seen in the 
background-subtracted PN and MOS-2 spectra 
shown in Fig.[5j The appearance of the dip may be 
due to unresolved strong emission lines. However, 
the wavelength of the dip coincides well with th e 
K-shell edges of oxygen (jVerner fc Yakovlevll995[ ). 
where our stellar wind models predict an absorp- 
tion maximum (see Fig. [2]). Therefore it is tempt- 
ing to attribute the observed dip to the oxygen 
K-shell absorption - to confirm this identification, 
data of higher quality will be required. Assum- 
ing the presence of K-shell of oxygen, the thermal 
plasma model fit to our low S/N data indicates 
presence of plasma with temperatures spanning 
from 1MK to 100 MK. 

4. On The Origin of X-ray Emission from 
WR 142 

The confident detection of weak hard X-rays 
from WR 142 prompts us to re-consider the pre- 
vious non-detections of X-ray emission from WR 
stars. The X-ray luminosity of WR 142 of Lx ~ 
7x 10 30 ergs -1 is comparabl e to present upper lim- 



its for the non-detections (lOskinoya et aD 2003; 



Oski noval l2005t iGosset et al.l 120051 ; I Skinner et al 
2006). Perhaps some WR stars that are yet un- 



detected in X-rays may be weak sources similar 
to WR142. This raises the question, what is the 
mechanism responsible for the generation of X- 
rays in WR stars? Usual suspects are wind shocks 
and magnetic fields. 

The strong-shock condition predicts the peak 
temperature in the shocked gas, Tx, as 



X 
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(1) 



where U is the velocity with which the gas rams 
into the shock (i.e., pre-shock velocity relative to 
the shock) , and ji is the mean molecular weight per 
particle in the post-shock gas. As the shock tem- 
peratures are very high, abundant species will be 
nearly entirely ionized. Therefore, fi depends only 
weakly on chemical compositions. For example, 



(i rs 4 f° r ionized H, while for fully ionized car- 
bon, fi « |. Thus, differences in wind molecular 
weight between O and WR stars cannot strongly 
affect the temperatures of the shocked plasma in 
stars of these types. 

Note that Eq. (JTJ) gives an upper limit to the 
temperature at the shock. Radiative cooling will 
restrict the hottest temperatures to a limited part 
of the total emission measure. Moreover, some 
fraction of the energy will be consumed by ion- 
ization processes. To roughly estimate the energy 
required to ionize WO wind material in the post- 
shock zone we consider oxygen and assume that 
the leading ion in the pre-shock material is O VI, 
while in the post-shock plasma it is O ix. A WO 
wind may contain 30% of oxygen. The ioniza- 
tion potentials are IP(0 Vl)« 0.14keV, IP(Q vn)» 
0.74 keV, and IP(Ovm)« 0.87keV dCoxl l2000h . 
The specific energy for ionization is e = EA^IPj, 
where TV; = Xi^Aima) -1 with X{ being the mass 
fraction of the element, A\ its atomic weight, and 
mu the atomic mass unit. Inserting these num- 
bers, the full ionization of oxygen requires e w 
3 x 10 13 ergg -1 . This is small compared to the ki- 
netic energy of the wind (5 x 10 15 ergg _1 for typi- 
cal 1000 km s -1 ). Thus the ionization is not a sig- 
nificant cooling process even in metal-rich winds. 

In th e time-dependent hyd rodynamic simula- 
tions bv lFeldmeier et aL ( 19971 ) the velocity jump 
U depends on the ratio between the period of 
the perturbations at the wind base, T c , and the 
flow time, Tfl ow = R^/v^. The former was esti- 
mated by the acoustic cut-off period, T c = aff -1 , 
where a is thermal speed and H is the pressure 
scale height. Assuming that the perturbation are 
seeded in the hydrostatic layers within the star, 



Tc/Tfio 



The thermal speed is lower 



in hydrogen-deficient WR atmospheres, while v esc 
can be larger than in O stars, as it is the case 
for our program star WR142. Therefore, the ve- 
locity jumps in WR wind shocks are expected to 
be rather smaller than larger, compared to O-star 
winds. 

These qualitative considerations indicate that 
if the same shock mechanism were in operation in 
WR and O star winds, the X-ray spectra from WR 
stars would be softer than those of O stars. How- 
ever, observations show the contrary - the handful 
of single WR stars with available spectral mea- 
surements all show X-ray emission harder than 
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typically found in O stars ( Skinner et al.|[2002af bl; 
Ignace et al 1 l2003h . This trend is confirmed by our 
new data on WR 142. 

Relatively hard X-ray emission can occur in a 
binary system; however, optical spectra of WR 142 
show no contamination from an early-type com- 
panion. A low-mass coronal companion to th is ini- 
tially very massive star is extremely unlikely (|Lucv 

20061). 

Babel fc Montmerld (|l997i) proposed that if a 



stellar wind is magnetically confined it can be 
strongly heated. The magnetic field locally dom- 
inates the bulk motions if the magnetic energy 
density exceeds the wind kinetic energy density, 
B 2 /p > pv 2 . At characteristic distance of 1 i?* 
from the photosphere, the velocity is v ~ 0.5voo ~ 
2000^8^, density p = M/A-Kv{r) 2 r « 1 x 
10- 10 gcm- 3 . Therefore, B(r = 2R*) > 7kG 
is required to control the WR 142 wind. For a 
dipole field, this implies that the magnetic field 
strength at the surface is larger than 50 kG. A 
field of such strength is not unrealistic, since the 
radius of WR 142 is ~ 50 times smaller than that 
of OB supergiants, where surface field s of 100 G 
have been observed i Bouret et al. I l2008h . One may 
speculate that even less strong magnetic fields may 
lead to hard X-rays due to some type of reconnec- 
tion and heating processes as invoked to explain 
the hi gh X-ray temperatures obs erved in O type 
stars (jWaldron fc Cassineln1l2007[) . 

It is difficult to test where in the wind the ob- 
served X-rays have been produced. Since no high- 
resolution X-ray spectrum of a single WR star 
exists, the f-i-r emission-line diagnostic that has 
been used for O-star winds has not been applied 
to WR stars. X-rays produced close to the stel- 
lar core must travel through the bulk of photo- 
absorbing wind. Tentative evidences for absorp- 
tion edges are foun d in the X-ray spectra of WR 1 
(jlgnace et al.ll2003f ) and WR142. 

In the wind of WR142, the radius at which 
the photo-absorption to X-rays is equal to unity is 
plotted in Fig. [2] against photon energy. The two 
shown models reproduce the stellar SED equally 
well. One of the models is thick to X-rays for 
1000's of and the observed X-rays would have 
to emerge from the far outer regions of the wind. 
Origin of X-rays far out in the wind is favored 
in the al ternative scen ario of X-ray emission from 
O stars ( Poll ock! 2007 ). However, a reduction in 



the mass-loss (a poorly constrained parameter) 
by a factor of only two and/or a higher effec- 
tive stellar temperature result in a higher degree 
of wind ionization. In this case a fraction of X- 
rays formed deep in the wind could freely escape. 
Win d clumping can further reduce wind attenua- 
tion (|Feldmeier et al.ll2003t ). 

There is suggestive evidence that WR 142 is 
a fast rotator, as is the case for the WN2 star 
WR2 (see Sect. [5]). The detection of strong and 
hard (log Lx » 10 32 ergs" 1 , T x > 10 MK) X-rays 



in WR2 (jSkinner et al.ll2008h likely indicates that 
a similar mechanism operates in these hot, com- 
pact, fast rotating stars. The lower X-ray lu- 
minosity of WR142 compared to WR2 can re- 
flect the higher absorption of X-rays in its denser 
wind. Fast rotation can lead to a slow equatorial 
wind with enhanced density, possibly forming an 
out-flowing disk, a nd a faster, thinner polar wind 
( Ignace et al. .1996). If the round shape of the line 
profiles found in both, WR 142 and WR 2, is due to 
rotational broadening, the inclination angle must 
be large. 

Summarizing, we have reported the detection of 
weak, but hard X-rays from WR 142 and speculate 
that their origin is connected with magnetic activ- 
ity of this far-evolved, compact WR star. The spe- 
cific shape of its optical line profiles signals that 
WR 142 may be rotating at nearly break-up ve- 
locity. In this case the spherical symmetry will 
be broken, potentially affecting the X-ray produc- 
tion and absorption. Observations of better qual- 
ity and progress in modeling are certainly needed 
to understand fully the high-energy processes in 
the winds of massive stars at late stages of their 
evolution. 

Based on observations obtained with XMM- 
Newton, an ESA science mission with instruments 
and contributions directly funded by ESA Mem- 
ber States and NASA. This research has made use 
of NASA's Astrophysics Data System Service and 
the SIMBAD database, operated at CDS, Stras- 
bourg, France. Special thanks to Robert Gru- 
endl for providing the IR measurements. The au- 
thors thank the referee M. De Becker for impor- 
tant comments that improved the paper. Fund- 
ing for this research has been provided by NASA 
grant NNX08AW84G (Y-HC and RI) and DLR 
grant 50 OR 0804 (LMO). 
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Fig. 1. — Spectral energy distribution of WR142. 
Boxes give observed 2MASS, Spitzer IRAC and 
MIPS magnitudes (labels). Observed optic and IR 
spectra are shown by thin blue lines. The red line 
is the continuum flux (plus lines in the UV and op- 
tical) from a PoWR model with logL/L Q = 5.35, 
T* = 160kK, M = lCT^-Mayr" 1 reddened with 
Eb-v = 1-73 mag. The model with logL/L© = 
5.65 and M = 10 _5 4 Af©yr _1 gives an equally 
good fit. A typical Hen line is shown in the inset 
to illustrate the effect of extreme rotational broad- 
ening - the observed line (blue) is compared to the 
synthetic line convolved with v sin i = 4000 km s _1 
(red). 
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Fig. 2. — Radius (in units of i?*) where the ra- 
dial optical depth becomes unity, as a function of 
wavelength in the X-ray range. The red curve is 
for the same model for which the spectral energy 
is shown in Fig.[TJ while the blue curve is for a 
model with half the mass-loss rate and otherwise 
similar parameters. 
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Fig. 4.— Thick red line: XMM-Newton MOS-2 
spectrum of WR 142 where background has not 
been subtracted. Thin blackline: MOS-2 spec- 
trum of a background region. The error bars give 
2a. 




Fig. 3.— Part of the merged XMM-Newton EPIC 
image (0.2-12.0 keV) with over-plotted contours. 
Image size is 5.7' x 4.1'. WR142 is marked by 
an arrow. The coordinates are equatorial (J2000). 
North up, east left. 
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Fig. 5.— XMM-Newton MOS-2 (red) and PN 
(blue) spectra of WR 142 with subtracted back- 
ground. The error bars give 2a 
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